        %% LoopSteps

        [Z_x,density_x]=DPR(Tpj,Ppj);
        Q_a=Q_cal(Z,Q_s,Tpj,Ppj);   %Nm3/s -> m3/s
        [lambda,a_coeff]=lambda_cal(Tpj,Ppj,d_x,Q_s,relative_den);
        %(Tpj,P1,lambda,a_coeff,s,Q,d,Z,relative_den)
        [P2,Ppj]=Colebrook(Tpj,P1,lambda,a_coeff,del_gc(lc_x),Q_s,d_x,Z_x,relative_den);
        %(T0,Tq,Pq,Pz,Qa,density,d)
        [T2,Tpj]=Temp_cal(Tground(1),T1,P1,P2,Q_a,relative_den,D_x);